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Abstract. We discuss the microwave energy losses in superconductors in the critical state. The field-induced 
variations of the surface resistance are determined, in the framework of the Coffey and Clem model, by 
taking into account the distribution of the vortex magnetic field inside the sample. It is shown that the 
effects of the critical state cannot generally be disregarded to account for the experimental data. Results 
obtained in bulk niobium at low temperatures are quantitatively justified. 

PACS. 74. 25. Ha Magnetic properties - 74.25.Nf Response to electromagnetic fields (nuclear magnetic 
resonance, surface impedance, etc.) - 74.60.Ge Flux pinning, flux creep, and flux-line lattice dynamics 



1 Introduction 



Investigation of fluxon dynamics in type-II superconduc- 
tors is of great interest for both fundamental and applica- 
tive aspects. From the basic point of view, it allows mea- 
suring the relative magnitude of elastic and viscous forces, 
which rule the motion regime of the fluxon lattice 
EE]. From the technological point of view, it allows de- 
termining the critical current, accounting for the energy 
losses and investigating the presence of irreversible phe- 
nomena, which are important factors for the implementa- 
tion of superconductor-based devices [S] . 

A suitable method to investigate the fluxon dynamics 
consists in determining the field-induced energy losses at 
microwave (mw) frequencies, by measuring the surface re- 
sistance, R s 0. In the absence of static magnetic fields, 
the variation with the temperature of the condensed-fluid 
density determines the temperature dependence of R s . On 
the other hand, the field dependence of R s in supercon- 
ductors in the mixed state is determined by the presence 
of fluxons, which bring about normal fluid in their cores, 
as well as the fluxon motion PJ2,3,4,6,7|H]. 

Coffey and Clem (CC) have elaborated a comprehen- 
sive theory for the electromagnetic response of type-II su- 
perconductors in the mixed state [H], taking into account 
flux creep, flux flow, and pinning effects in the framework 
of the two-fluid model. The theory has been developed 
with two basic assumptions: i) inter- vortex spacing much 
less than the field penetration depth; ii) uniform vortex 
distribution in the sample. With these assumptions, the 
local vortex magnetic field, B(r), averaged over distances 
larger than several inter- vortex spacings, is spatially uni- 
form. As pointed out by the authors [Sj, the first assump- 
tion is valid when the external field, Hq, is greater than 



2-ffei; whereas, the second one does not take into account 
the magnetic history of the sample. 

It is well known that the H - T phase diagram of 
type-II superconductors is characterized by the presence 
of the irreversibility line, Hi rr (T), below which the mag- 
netic properties of the superconducting sample become ir- 
reversible. The application of a DC magnetic field smaller 
than Hi rr (T) develops a critical state of the vortex lat- 
tice l^lllOj and the assumption that B is uniform over the 
sample is not longer valid. Although the critical state is a 
metastable energy state, the relaxation times toward the 
thermal equilibrium state can be much longer than the 
time during which the measurements are performed 
In this case, the effects of the critical state can be revealed 
in all the superconducting properties involving the pres- 
ence of fluxons in the sample. 

When the sample is in the critical state, the CC model 
does not account for the experimental results; for instance, 
no magnetic hysteresis in the R s (Ho) curves is expected. 
Though different authors have justified the hysteretic be- 
havior of the R s (Hq) curves by considering critical-state 
effects in the fluxon lattice |12U1."{| . the field dependence 
of the surface resistance for non-uniform flux profile has 
never been quantitatively investigated. 

In this paper, we show that, when the superconduct- 
ing sample is in the critical state, the CC theory has to 
be generalized by properly taking into account the flux 
distribution inside the sample. We suggest a way to con- 
sidering the critical-state effects. We will show that the 
R S (H ) curve strongly depends on the specific profile of 
the magnetic induction B, determined by the field depen- 
dence of the critical current density J C {B). We will report 
experimental results of the field-induced variations of R s 
in a niobium bulk sample in the critical state; the results 
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are quite well accounted for by considering a specific B 
profile inside the sample. 



2 The model 

In the London local limit, the surface impedance is pro- 
portional to the complex penetration depth A of the em 
field. In particular, 

R 8 = -fi ujIm[X(cj,B,T)]. (1) 

In the CC model, X(u>,B,T) is given by 



r, p m_ X^B,T) + (i/2)S^,B,T) 



with 



X{B,T) 



Ao 



6 nf {u,B,T) 



y/[l-(T/T c ni-B/B c2 (T)Y 



So 

^/1-{1-(T/T c ni-B/B C2 (T)} 



(3) 



(4) 



where Ao is the London penetration depth at T — and 
5q is the normal- fluid skin depth at T = T c . 
5 V is the effective complex penetration depth arising from 
the vortex motion and depends on the relative magnitude 
of the viscous-drag and restoring-pinning forces. An im- 
portant parameter, which determines the regime of the 
fluxon motion, is the so-called depinning frequency, lu c , de- 
fined by the ratio between the restoring-force and viscous- 
drag coefficients. When the frequency of the em wave, u>, 
is much larger than oj c the fluxon motion is ruled by the 
viscous-drag force and the induced em current makes flux- 
ons move in the flux-flow regime, even for Hq < Hi rr (T). 
On the contrary, for lo -C u) c the motion of fluxons is ruled 
by the restoring-pinning force. In the following, we will de- 
vote the attention to the case in which the vortex lattice 
moves in the flux-flow regime, where the field-induced en- 
ergy losses are noticeable. In this case, considering the 
expression of the viscous coefficient proposed by Bardeen 
and Stephen 14 , it results 6% = 5%B/B c2 (T). The real 
part of A defines the AC penetration depth, X ac - 

Since the mw losses depend on the local magnetic field 
B(r), when it is not uniform over the sample, different re- 
gions of the sample contribute to the energy losses differ- 
ently; this should occur when the sample is in the critical 
state. Generally, the critical state develops at tempera- 
tures smaller enough than T c , where the pinning effects 
are significant; in this case, the energy losses are mainly 
related to the motion of fluxons induced by the mw cur- 
rent. So, a discriminating criterium, to determine in what 
extent the non-uniform B profile affects the energy losses, 
consists in evaluating the variation of the DC magnetic in- 
duction in the region where fluxons feel the Lorentz force. 



Fig. n describes the distribution of the fields and the 
mw current in the two geometries: H U \\H$ (a); H^-LHo 
(b). For the sake of clearness, we limit the analysis to a 
slab geometry (sample width w and sample height d, with 
X ac <C w, d). We suppose that the static magnetic field, 
Hg, is greater than H c i and that the sample is in a critical 
state a la Bean, i.e. characterized by a field-independent J c 
[5|; however, the analysis applies to any J C (B) dependence. 
When the mw magnetic field is parallel to the z-axis, the 
mw current penetrates in the surface layers of the sample, 
of width X ac , in the x - y plane. If H u \\Ho, only the fluxons 
in these layers experience the Lorentz force, due to the mw 
current, and the fluxon lattice undergoes a compressional 
motion. In Fig. 1 (a) we analyze the regions indicated by 
the shadowed area. In this case, if the variation of B in 
these regions is negligible, i.e. J c X ac <C Hq, the vortex 
magnetic field can be considered as uniform. A different 
situation occurs when H^J-Hq, as shown in Fig. ^(b); in 
this case, all the vortices present in the sample are involved 
in the motion, the fluxon lattice undergoes a tilt motion 
[7J, X ac determines the part of the flux line in which the 
Lorentz force acts; the approximation of B uniform would 
be proper only if the sample width is such that J c w <C Hq . 
So, for bulk samples, it is expected that the specific profile 
of B strongly affects the R S (H ) curve; indeed, different 
regions of the sample contribute in different extent to the 
mw energy losses. In the following, we will refer to this field 
geometry, where the effects of the non-uniform profile of 
B most affect the R s (Ho) curve. 

Brandt [7] has shown that the compressional and tilt 
motion of the fluxon lattice can be described by the same 
formalism and the same AC penetration depth results be- 
cause the compression and tilt moduli are approximately 
equal. So, the results obtained from the CC model are 
valid even in the case of tilt motion provided that B is 
uniform over the region where the Lorentz force is active. 

When the DC magnetic induction is not uniform, one 
can subdivide the sample surface in different regions, in 
each of them B(r) can be considered uniform. Each part 
of the surface is characterized by a different R s value, due 
to the local magnetic induction, and the energy losses of 
the whole sample are determined by the surface-resistance 
contributions of each region. The measured surface resis- 
tance is an average over the whole sample: 



fl s = ~ J R s (\B(r)\)dS, 



(5) 



where £ is the sample surface, S is its area and r identifies 
the surface element; the main contribution comes from the 
sample regions where H x J w ^ 0. 

Usually, in order to disregard the geometrical factor of 
the sample, the reported experimental results of the sur- 
face resistance are normalized to the normal-state value 
at T = T c , R n . In the flux-flow regime, in order to cal- 
culate R s /R n , by Eqs. (1-4), it is necessary to know the 
ratio Xq/Sq and the upper critical field. Moreover, to take 
into account the critical-state effects, by Eq. it is also 
necessary to know the B profile inside the sample, deter- 
mined by J C (B). We would remark that the analysis can 
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Fig. 1. Distribution of the fields and current in two differ- 
ent geometries: fl^||J3o (a); H^-LHo (b); the shadowed area 
indicates the regions of the sample interested by the vortex 
motion. 



be generalized, in the different motion regimes, by intro- 
ducing the field dependence of the depinning frequency 
in 5 v (u>, B,T); however, this introduces further fitting pa- 
rameters. 

Fig- HI shows R s /R n expected in the flux-flow regime 
as a function of the reduced field, Hq/H C 2, in three dif- 
ferent cases. Curve (a) is the one expected from the CC 
model (J c = 0). Curves (6) and (c) have been obtained 
with J c = J c o, independent of B, and two different values 
of the full penetration field H* [§]. All the curves have 
been obtained with T = T c /4 and Xq/Sq = fCP 2 ; further- 
more, for simplicity, it was supposed H c \ = 0. ft is evident 
that, taking into account the distribution of B, different 
field dependencies of R s arise, fn particular, a comparison 
between curves (b) and (c) shows that the greater H* , the 
lower the value of R s ; this occurs because the induction 
field of the whole sample decreases on increasing H* . 

Another important characteristic of the R s (Hq) curves 
is the change of concavity occurring at Hq = H*. When 
the fluxon lattice moves in the flux-flow regime and B is 
uniform, it is expected R s cx \/B, with B « noHo, which 
leads to a negative concavity of the R S (H ) curve. For 
the general case of spatially-dependent flux density, the 
shape of the R s (Hq) curve is determined by the external- 
field dependence of B. fn particular, when in the whole 
sample the local magnetic induction linearly depends on 
the external field a negative concavity of the R s (Ho) curve 
is expected. This occurs in the thermal equilibrium state 



for Hq > 2H c i, as well as in the critical state a la Bean for 
Hq > H* [Sj. The positive concavity of the curves (b) and 
(c) of Fig. 2, obtained for Hq < H*, can be qualitatively 
understood considering that, on increasing the external 
field from H c \ up to H* , more and more regions contribute 
to the mw losses. 




o.o 



0.0 0.1 0.2 0.3 0.4 0.5 

H IH c2 

Fig. 2. Normalized R s (Ho) curves, expected in the flux-flow 
regime, in the different cases: (a) B uniform; (b) critical state 
a la Bean with H* = # c2 /40 and (c) H* = # c2 /4. T = T c /4, 
Ao/5 = 1(T 2 . 




Fig. 3. Expected R s (Ho) curves, obtained using for the field 
dependence of the critical current density the laws shown in 
the inset. T = T c /4, X /S = 10" 2 . 

When the field dependence of the critical current is 
taken into account, the shape of the R s (Hq) curve strongly 
depends on the B profile due to the specific J C {B) behav- 
ior. The field dependence of J c has been widely discussed 
in the literature both theoretically [fOlll 5llf 6j and exper- 
imentally |17U18U19j : different behaviors have been sug- 
gested. In Fig. we report the expected R S (H ) curves, 
obtained considering different J C (B) laws, along with the 
curve expected from the CC theory (a); also in this case, 
we have assumed H cl = 0. Curve (b) has been obtained 
using J c = Jo exp(— B/Tq), curve (c) with J c = Jq — a B 
and curve (d) with J c — J ; the inset shows the corre- 
sponding J C (B) laws used. As one can see, the features 
of the R s (Hq) curve strongly depend on the profile of B 



4 



M. Bonura et al.: Critical-state effects on microwave losses 



inside the sample. In particular, looking at curve (&), it is 
evident that when the critical current verges on zero the 
R s (Hq) curve approaches the one obtained from the CC 
theory, as expected. This finding is expected, in any case, 
for magnetic fields close to H C 2(T), where the pinning be- 
comes ineffective and, consequently, J c « 0. 

It is worth noting that we have investigated the effects 
of the B profile determined by the field dependence of the 
critical current density on the R s (Hq) curve, neglecting 
possible variation of the profile due to the finite dimen- 
sions of the superconducting sample. As it is well known, 
the sheet current, induced near the sample edge's, can af- 
fect the B profile in the sample (20j • These effects are par- 
ticularly enhanced for thin strips when the DC magnetic 
field is parallel to the smallest dimension; in this 
proper B profile has to be considered. 

3 Comparison with experiments 

In order to study the peculiarities of the R s (Hq) curve in 
superconductors in the critical state, we have measured 
the field-induced variations of R s in bulk Nb at low tem- 
peratures, where the pinning effects are significant and a 
magnetic field smaller enough than H c2 (T) develops the 
critical state of the fluxon lattice. The sample, of dimen- 
sions 3.3 x 2.3 x 1.3 mm 3 , has been cut from a Tokyo 
Denkai batch with RRR = 300 and undergoes a super- 
conducting transition at T c w 9.2 K. We have studied the 
Nb because the depinning frequency |23| is small enough 
to suppose that the mw current induces the fluxon lattice 
moving in the flux-flow regime; furthermore, our experi- 
mental apparatus allows reaching DC magnetic fields of 
the order of H c2 . 

The mw surface resistance has been measured using 
the cavity-perturbation technique |21| . A copper cavity, of 
cylindrical shape with golden-plated walls, is tuned in the 
TEon mode, resonating at 9.6 GHz. The cavity is placed 
between the poles of an electromagnet, which generates 
magnetic fields up to ~ 1 T. The sample is located in 
the center of the cavity where the mw magnetic field is 
maximum. The field geometry is that shown in Fig. 1 (b) 
and, therefore, the mw current induces a tilt motion of 
the whole fluxon lattice. The R s values are determined 
measuring the variation of the quality factor of the cavity, 
induced by the sample, by an /ip-8719D Network Analyzer. 

In order to check the presence of the critical state in 
the sample, we have measured the field-induced variations 
of R s by cycling the DC magnetic field from zero to the 
maximum value available in our experimental apparatus 
and back, at different values of the temperatures. Up to 
2 K below T c , R s exhibits a magnetic hysteresis that dis- 
appears for Hq higher than a certain value, depending on 
the temperature. In particular, in the range of tempera- 
tures T ss 2 -f- 4 K the hysteretic behavior is detectable 
up to Hq ~ 0.4 T; at higher temperatures, the field range 
in which the hysteresis is present shrinks; eventually, for 
T > 7 K, the R s (Ho) curve becomes reversible in the 
whole range of fields investigated. These findings confirm 
that, at temperatures smaller enough than T c , the external 



field develops a critical state in the fluxon lattice. Further- 
more, in the field range in which the hysteresis is observed, 
R s does not exhibit any detectable time evolution, from 
the instant in which the field value has been set, ensur- 
ing that the critical state does not relax during the time 
in which the measurements are performed. The hysteretic 
behavior of microwave surface resistance will be discusses 
in a forthcoming paper; in this paper, we are interested to 
investigate the critical-state effects on the R s (Ho) curve 
at increasing magnetic field, and to test the method we 
suggest for taking into account the not uniform B distri- 
bution in the sample. 

Fig. 4 shows the field-induced variations of R s , ob- 
tained in the Nb sample at the two temperatures T — 
4.2 K (a) and T = 2.2 K (b). In the figure, AR S (H ,T) = 
R s (Hq , T) — R res , where R res is the residual mw surface 
resistance at T = 2.2 K and Hq = 0; the data are nor- 
malized to the maximum variation, AR™ ax = R n — R res . 
R s does not show any variation as long as Hq reaches a 
certain value that can be identified as the first penetration 
field, H p . From the data, we obtain 7J p (4.2 K) w 52 mT 
and H p (2.2 K) w 57 mT; these values are smaller than 
those reported in the literature for the lower critical field 
in Nb superconductor because of the demagnetization ef- 
fects. The dashed and pointed lines are the curves ex- 
pected from the CC model; the continuous lines are those 
obtained in the framework of our model; all the expected 
curves are plotted for Hq > H p . Since the depinning fre- 
quency reported in the literature for Nb are smaller than 
100 MHz [23! > we nave assumed that fluxons move in the 
flux-flow regime. The essential parameters to calculate the 
expected curves are: Xo/Sq, iJ c2 (T) and the B profile in- 
side the sample. According to the results reported in the 
literature we have used Ao/<5o = 3 x 10~ 2 ; how- 

ever, for T < T c /2 the expected results are little sensitive 
to variations of this parameter. _ff C 2(4.2 K) has been de- 
termined experimentally and it results LIT; the value of 
the upper critical field at T = 2.2 K has been used as 
fitting parameter and we have found H c2 {2.2 K) = 1.8 T. 
To take into account the presence of the critical state, 
we have used a field dependence of J c similar to that re- 
ported by Cline et al. for Nb crystals; in particular, it 
has been used a linear field dependence of J c at low fields, 
followed by an exponential decrease for Hq > H 1 . The 
value of H' has been considered as fitting parameter; the 
best-fit curves (continuous lines of Fig. 4) have been ob- 
tained using H'(4.2 K) = 0.35 T and H'{2.2 K) = 0.5 T. 
The B profiles, obtained on increasing Hq, which come 
out from the used J C (B) laws, are shown in the insets. 
The full penetration fields calculated from the B profiles 
are H*(4.2 K) = 0.41 T and H*(2.2 K) = 0.56 T. 

As one can see, in the whole range of Hq investigated 
the experimental results are quite well accounted for by 
considering the B distribution inside the sample. At low 
fields, the non-uniform B profile affects to a detectable 
extent the energy losses, giving rise to the positive con- 
cavity of the R s (Hq) curve. On increasing Hq, the slope of 
the field profile decreases; the effects of the non-uniform 
B distribution become less and less important; eventu- 
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Ho"o(T) 

Fig. 4. Normalized field-induced variations of R s at: T — 
4.2 K (a) and T = 2.2 K (b). Symbols are the experimen- 
tal results. Dashed and pointed lines have been obtained from 
the CC model using \ /5 = 3 x 10" 2 ; H c2 (2.2 K) = 1.8 T 
and i/ C 2(4.2 K) = 1.1 T. The continuous lines are the best-fit 
curves of the data, obtained using the same values of Ao/<5o 
and H C 2(T) and the B profiles shown in the insets. 

ally, when B becomes roughly uniform over the sample, 
the observed behavior of R s (Hq) is consistent with that 
expected from the CC model. 

It has been experimentally observed and theoret- 
ically justified |27| . that, in samples of finite dimensions, 
the application of an AC magnetic field normal to the 
DC field gives rise to a "shaking" of the fluxon lattice, 
inducing relaxation toward the uniform distribution. The 
process is particularly relevant when the amplitude of the 
AC field is of the order of the full penetration field. In 
mw measurements with the cavity perturbation technique, 
the amplitude of the mw magnetic field is, usually, much 
smaller than 1 Gauss, so the process can play a role only 
in very thin samples or in proximity of the irreversibility 
line, where the critical current is very small. Of course, if 
this process occurs, the mw field assists fluxon distribut- 
ing uniformly in the sample and, consequently, a variation 
of the R S (H ) curve is expected. 

Another effect we have neglected is the deformation of 
the B profile near the edges of the sample [20]; however, 
the measured R s value is due to the average response of 
the whole sample, and, therefore, our measurements do 
not allow estimating in what extent this effect influences 
the experimental curves. On the other hand, the results 
reported in Fig. 4 show that the B profiles shown in the 
insets are suitable to well account for the experimental 
data. 



4 Conclusion 

In summary, we have studied the field-induced energy 
losses in superconductors in the critical state. We have 
shown that the distribution of the local vortex magnetic 
field in the sample can strongly affect the field dependence 
of the mw surface resistance. The analysis has been car- 
ried out for the field geometry in which the mw current 



induces a tilt motion of the fluxon lattice; indeed, in this 
case, the non-uniform flux distribution most affects the 
R s (Hq) curve. The expected curves have been obtained 
supposing that the fluxons move in the flux-flow regime, 
where the field-induced energy losses are relevant even for 
applied fields much smaller than H C 2] on the other hand, 
the analysis can be easily extended to a more general case, 
provided that the field dependence of the depinning fre- 
quency is known. We have highlighted that the effects of 
the critical state on the mw energy losses cannot gener- 
ally be disregarded to account for the experimental data. 
Experimental results of field-induced variations of R s in a 
Nb bulk sample at low temperatures have been justified 
quite well in the framework of our model. 

The authors are very glad to thank G. Lapis and G. 
Napoli for technical assistance. 
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